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ABSTRACT 


Using a 2.5-dimensional MHD simulation, we investigate the role played by the inner coronal null point in the 
formation and evolution of solar quiescent prominences. The flux rope is characterized by its magnetic fluxes, the 
toroidal magnetic flux ®, and the poloidal flux &5,. It is found that for a given ®,, the catastrophe does not occur 
in the flux rope system until G5, increases to a critical point. Moreover, the magnetic flux of the null point is the 
maximum value of the magnetic flux in the quadrupole background magnetic field, and represented by wy. The 
results show that the bigger yy usually corresponds to the smaller catastrophic point, the lower magnetic energy 
of the flux rope system, and the lesser magnetic energy inside the flux rope. Our results confirm that catastrophic 
disruption of the prominence occurs more easily when there is a bigger wy. However, yy has little influence on 
the maximum speed of the coronal mass ejections (CMEs) with an erupted prominence. Thus we argue that a 
topological configuration with the inner coronal null point is a necessary structure for the formation and evolution 
of solar quiescent prominences. In conclusion, it is easier for the prominences to form and to erupt as a core part 
of the CMEs in the magnetic structure with a greater wy. 


doi:10.1088/0004-637X/800/1/43 
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1. INTRODUCTION 


The formation and evolution of solar quiescent prominences 
(or filaments) is an important unsolved problem in solar physics. 
At present, it is widely accepted that the magnetic field plays 
an important role in the filament eruption process. The original 
model for solar prominences relied on dips in the magnetic 
to hold the cool plasma against gravity (e.g., Kippenhahn & 
Schlüter 1957; Hyder 1967; Kuperus & Raadu 1974; Heyvaerts 
et al. 1977; Spicer 1977; Milne et al. 1979; Anzer & Priest 
1985; Zirker 1989; Antiochos et al. 1994; Forbes et al. 1994; 
Low & Hundhausen 1995; Tandberg-Hanssen 1995; Aulanier & 
Demoulin 1998; Martin 1998; Heinzel & Anzer 1999; Labrosse 
et al. 2010; Mackay et al. 2010; Low et al. 2012). This classic 
model has been extended in three ways to account for energy 
storage and release. The first considers magnetic reconnection 
in the extended current sheet and on the tearing mode instability 
by magnetic shear in the dynamics of prominences (e.g., Foukal 
& Behr 1995; DeVore & Antiochos 2000; Petrie & Low 
2005; Welsch et al. 2005; Luna et al. 2012; Xia et al. 2012). 
The second uses the cancelation of converging magnetic flux 
(e.g., Galsgaard & Longbottom 1999; Martens & Zwaan 2001; 
Litvinenko & Wheatland 2005; López et al. 2006; Mackay & 
van Ballegooijen 2006a, 2006b; Litvinenko et al. 2007). The 
third relies on the emergence of a magnetic flux rope (e.g., 
Low & Zhang 2002; Fong et al. 2002; Aulanier & Demoulin 
2003; Low & Petrie 2005; Okamoto et al. 2009; Magara 2012; 
Hillier & van Ballegooijen 2013; Zhang 2013; Longcope & 
Forbes 2014; Kliem et al. 2014). In this work, we discuss the 
dynamical evolution of a flux rope system in a complicated 
background field. 

Observations show that the prominences are usually sus- 
pended in the inner corona within a magnetic complex consisting 
of regions (e.g., Newkirk 1967; Saito & Hyder 1968; Bohlin 
1970; Waldmeier 1970; Pneuman 1972; Saito & Tandberg- 
Hanssen 1973; Webb et al. 1976; Schmahl & Hildner 1977; 
Serio et al. 1978; Dryer 1982; MacQueen et al. 1983; Illing 


1984; Athay & Illing 1986; Iling & Hundhausen 1986; 
Gopalswamy & Kundu 1990; Raju et al. 1993; Gopalswamy 
et al. 1997; Wiik et al. 1994; Low 2001; Marqué 2004; 
Koutchmy et al. 2004; Gibson et al. 2006, 2010; Wang et al. 
2006; Zhou et al. 2006a, 2006b; Heinzel et al. 2008; Zhang 
& L1 2009; Habbal et al. 2010; Su et al. 2010; Li et al. 2011; 
Régnier et al. 2011; Berger et al. 2012; Vourlidas et al. 2013; 
Parenti 2014; Yang et al. 2014). Introducing the quadrupole 
field to flux rope catastrophe model, we found that there are two 
current sheets in the flux rope system, one is a vertical current 
sheet connected to the solar surface and the flux rope, and the 
other is a transversal sheet evolved from the null point, which 
was originally located in the inner corona (Zhang et al. 2005; 
Zhang & Wang 2007). Furthermore, the recent paper presented 
by Zhang (2013) shows that after the first catastrophe the qui- 
escent prominences are supported by the vertical and transverse 
current sheets in the inner corona, and after the second catas- 
trophe the prominence would erupt as a part of coronal mass 
ejections (CMEs). The magnetic null point is an important part 
of the magnetic field structure, and it is essential to the forma- 
tion of the current sheet. Therefore, the purpose of this paper is 
to study the role of the inner corona null point in the formation 
and eruption of quiescent prominences using a magnetic flux 
rope catastrophe model. The model is introduced in Section 2, 
and the numerical results are discussed in Section 3. Finally, in 
Section 4, we discuss several conclusions. 


2. NUMERICAL MODEL 


Using time-dependent resistive MHD simulations, we study 
the evolution of solar quiescent prominences. For 2.5- 
dimensional MHD problems in spherical coordinates (r, 0, o), 
a magnetic flux function y(t, r, O) was introduced, and it is 
related to the magnetic field by 


y. i 
B= = B B, =B 1 
V x (Au + Q> p gP, ( ) 
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Figure 1. Magnetic configurations of the flux rope system corresponding to two cases: for case A, a = 1.0. P, = 0.36 and two separate values of By: (a) 0.0621 right 


before and (b) 0.0622 right after the first catastrophic point. For case B, œ = 1.01. ®, = 0.36 and two separate values of My: (c) 0.0641 right before and (d) 0.0642 
right after the first catastrophic point. The outline of the electric current sheets is depicted by thick solid curves. 


where B, is the azimuthal component of the magnetic field. 
Then the corresponding resistive MHD equations are cast in the 
following form: 

p 


p ee (2) 
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y (= 1.05) is the polytropic index, 7 is the resistivity, and j is the 
current density. The solution domain is taken to be 1 < r < 30, 
O < 0 x 7/2, discretized into 130 x 90 grid points. The 
MHD equations are solved with the multi-step implicit scheme 
developed by Hu (1989). 

The initial quiet corona is assumed to be isothermal and static 
with Ty = 1.4 x 106 K and pọ = 6.68 x 107? kg m^? at the 
coronal base. The ratio of gas pressure to magnetic pressure is 
B = 0.03, and then the unit of y is Wo = (2upoRTo R4/p)'/? = 
5.50x10!^ Wb, the unit of magnetic field intensity is Bo = 
Vo/ R2 = 1.14 x107? T, va = Bo/(upo)!? = 1242 kms^! for 
velocity, r4 = R./v4 = 560s for time, the unit of magnetic 
energy is Eo = BG R;/p = 3.47 x10% J, and jo = Bo/(u Ro) 
= 1.31x1079A m^? for current density. 

The initial background magnetic field is a quadrupole field 
(Antiochos et al. 1999; Hu 2004; Zhang et al. 2005; Zhang 
2013), and the related magnetic flux function normalized by yo 
is given by 


vr, 0) = 


sin?6 (3 +5cos 20) sin? 0 
+ (9) 
" 


2r3 , 


where r is in units of R,, and o is the control parameter of 
magnetic flux. Moreover, the magnetic flux of the null point 
is the maximum value of magnetic flux in the quadrupole 
background magnetic field, and is represented by wy. During 
a magnetic flux rope emerging from the bottom of the corona 
into the central arcade, the null point evolves into the transverse 
current sheet. 
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Figure 2. For case A: heights of the coronal magnetic flux rope vs. time (panels (a) and (b)), velocity of the flux rope axis vs. time (panels (c) and (d)), and magnetic 


energy of the system and of the flux rope vs. time (panels (e) and (f)). 


Furthermore, we calculate the magnetic energy E of the force- 
free field, normalized by 47 Eo, using 


] 20 z/2 
E= >| ar f B’r? sin6d6 
2 Ji 0 


30? z/2 
p Bi 


7 (B? — Bj) „3 Sin 0d9, 


(10) 


where the first term on the right-hand side is the magnetic energy 
in the computational domain (1 < r < 30; 0 < 0 < mz; 
0 < @ < 2x) and the second is that outside of the domain, 
having been transformed into a surface integral over the top (Hu 
2004). In addition, magnetic energy inside the flux rope, E,, is 
obtained by an integral representation similar to the first term 
on the right-hand side of the Equation (10), but with the integral 
domain limited to the interior of the rope (y < 0). 


Zhang (2013) shows that the formation of solar quiescent 
prominences may well be an ideal MHD process, and the 
prominences are supported by the vertical current sheet and 
located in cavities below the curved transverse current sheet 
in the inner corona after the first catastrophe. However, in the 
process of prominence eruption, magnetic reconnection occurs 
in the vertical and transverse current sheets after the second 
catastrophe. Moreover, we introduce a critical current density 
for the two current sheets, the vertical one is j, — 20 and the 
transverse one is j; = 5 (see Zhang & Wang 2007). Besides, 
we choose the critical length of the vertical current sheet, 
le = 0.5 Ro, which is large enough to satisfy the tearing mode 
instability (Furth et al. 1963). Therefore, the current density 
nearby the transverse current sheet exceeds j;, or that nearby the 
vertical current sheet exceeds j, and the length of the vertical 
current sheet also exceeds l+; n is set to 0.01, and 77 is set to zero 
elsewhere (Zhang et al. 2006). In the present study, the flux rope 
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Figure 3. Magnetic configuration in black solid curves and temperature distribution in color at several separate times for case A, in which just before and after the 
first catastrophe (panels (a) and (b)), just before and after the second catastrophe (panels (c) and (d)), reconnection occurs in the vertical and transverse current sheets 
(panel (e)), and then the velocity of the CMEs reaches the maximum speed (panel (f)). The cool blue and hot red correspond to 1.58 x 10° and 2.51 x 10°K, 


respectively. 


is characterized by the toroidal magnetic flux ®, and poloidal 
flux 4. First, in the ideal MHD simulation, it is found that for a 
given ® ,, the catastrophe does not occur in the flux rope system 
until ®, increases to a critical point. ®, is the catastrophe point 
for ®,. Then in the time-dependent resistive MHD simulation, 
the flux rope represented by the ®, and the P, has emerged 
from the bottom of corona into the central arcade. The evolution 
of the flux rope system shows that after the first catastrophe 
the prominences are formed in the inner corona, and after a 
period of time the second catastrophe occurs, meanwhile, the 
prominences erupt as a main part of CMEs that resulted in the 
magnetic reconnection taking place in the two current sheets. 


3. SIMULATION RESULTS 


In our simulation, two parameters are used to describe 
the geometrical characteristics of the flux rope system in 
equilibrium. One is h,, the length of the vertical current sheet, 
the other one is h,, the height of the quiescent prominences 
(HQP). Moreover, the magnetic properties of the flux rope 
are represented by ®,, the toroidal magnetic flux, and ®,, 
the poloidal magnetic flux. In order to investigate the role of 
the inner coronal null point in the formation and evolution of 
quiescent prominences, two different o are chosen: 1.0 and 1.01. 


Table 1 
Magnetic Energies of the Flux Rope System and Magnetic Flux of the 
Maximum in the Magnetic Configuration, wy, for Each œ 


a VN Ep 
1.0 0.3849 1.476 
1.01 0.3830 1.499 


Note. See the text for details. 


Then, we find two cases: case A and case B, respectively. When 
a increases, Wy decreases, and the magnetic energy of the flux 
rope system Eo increases (see Table 1). 


3.1. Formation of Solar Quiescent Prominences 


In the quadrupole magnetic field structure, quiescent promi- 
nences may be formed by the first catastrophe (Zhang 2013). In 
this simulation, with the same emerging strength of magnetic 
flux, P, = 0.36, we obtain two different catastrophic points, 
P, = 0.0622 and 0.0642, for which œ = 1.0 and 1.01, respec- 
tively. For case A, at P, = 0.0621, the flux rope sticks to the 
solar surface, and a curved transverse current sheet is formed 
above the flux rope (Figure 1(a)). Then at P, = 0.0622, the flux 
rope breaks away from the solar surface, and suspends stably 
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Figure 4. For case B: heights of the coronal magnetic flux rope vs. time (panels (a) and (b)), velocity of the flux rope axis vs. time (panels (c) and (d)), and magnetic 


energy of the system and of the flux rope vs. time (panels (e) and (f)). 


in the inner corona, for which the flux rope is located above 
the short vertical current sheet and below the transverse cur- 
rent sheet (Figure 1(b)). Similarly, for case B, at P, = 0.0641, 
there has been no catastrophe in the system (Figure 1(c)), and 
right at P, = 0.0642, the first catastrophe occurs (Figure 1(d)). 
Furthermore, comparing case A and case B, we find that the 
bigger Yy corresponds to the smaller catastrophic point, the 
lower magnetic energy of the flux system, the lesser magnetic 
energy inside the flux rope, and the similar HQP (see Tables 2 
and 3). 


3.2. Evolution of Solar Quiescent Prominences 


In the simulation, the magnetic flux rope represented by 
(®,, ®,) has emerged from the bottom of the corona into 
the central arcade, and we obtain the results of the evolution 
of flux rope system. With the purpose of studying the effect 


Table 2 
Geometrical Characteristics and Magnetic Energies of the Flux Rope System 
Just Before and After the First Catastrophe for Case A, for which o = 1.0 


$, ®, he HQP E E, 
(Rs) (Rs) 

0.36 0.0621 0.0 0.1656 1.709 0.1491 

0.36 0.0622 0.02209 0.1780 1.713 0.1497 


Note. See the text for details. 


of œ on prominence evolution, (®, = 0.36, P, = 0.0622) 
and (, = 0.36, P, = 0.0642) are chosen for case A and 
case B, respectively. Figure 2(a) shows that the flux rope has 
emerged into the arcade, and then the prominences are formed 
by the first catastrophe after r4 = 10. Both HQP and h, rise 
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Figure 5. Magnetic configuration in black solid curves and temperature distribution in color at several separate times for case B, in which just before and after the first 
catastrophe (panels (a) and (b)), just before and after the second catastrophe (panels (c) and (d)), reconnection occurs simultaneously in the vertical and transverse 
current sheets (panel (e)), and then the velocity of the CMEs reaches the maximum speed (panel (f)). The cool blue and hot red correspond to 1.58 x 106 and 


2.51 x 10° K, respectively. 


Table 3 
Geometrical Characteristics and Magnetic Energies of the Flux Rope System 
Just Before and After the First Catastrophe for Case B, for which a = 1.01 


9, o, he HQP E E, 
(Rs) (Rs) 

0.36 0.0641 0.0 0.1681 1.739 0.1577 

0.36 0.0642 0.02171 0.1791 1.745 0.1575 


Note. See the text for details. 


slowly until the second catastrophe occurs. After c4 = 29, HQP 
increases from 0.2345 R, to 0.6252R,, and A, increases from 
0.05422 R, to 0.3881 R,, respectively (Figure 2(b)). Because of 
the prominence eruption, HQP increases quickly. Figure 2(c) 
shows that the velocity of the flux rope was within the range of 
0.5 kms^! and 1.5kms~! before 20t, except that the velocity 
increased rapidly and then fell quickly before 5r4. Meanwhile, 
the velocity has no obvious change even though the first 
catastrophe occurs. Figure 2(d) shows that the prominences 
have an oscillation process before the second catastrophe, and 
then the prominences erupt to reach the maximum speed of 
1055 km s^! at 3274. In the beginning of the emergence of the 
flux rope, the magnetic energy of the flux rope system increases 
from 1.476 to 1.769 in a very short period of time. However, 
soon afterward, the energy decreases to 1.704, and then it has 


no change even though the first catastrophe occurs. Until the 
second catastrophe happens, it significantly decreases to 1.582 at 
which the prominences reach the maximum speed (Figure 2(e)). 
Furthermore, the variation curve of magnetic energy inside the 
flux rope is similar to that of the flux rope system (Figure 2(f)). 

Figure 3 describes the evolution of magnetic flux rope system 
at several separate times for case A. Before t = 10t,, the 
catastrophe does not occur for which the flux rope below 
the undeveloped transverse current sheet attaches to the solar 
surface (Figure 3(a)). At about t = 1174, the first catastrophe 
occurs, and the length of the vertical current sheet below the 
rope is 0.02007 R, (Figure 3(b)). The quiescent prominences 
have levitated in the inner corona until t = 29t, (Figure 3(c)). 
Then, at about t = 3074, the second catastrophe occurs and the 
prominences erupt, namely, the CMEs onset. However, at this 
time, magnetic reconnection does not occur in the elongated 
vertical current sheet and the curved transverse current yet 
(Figure 3(d)) Obviously, Figure 3(e) shows that magnetic 
reconnection has just happened in the transverse current sheet 
while it has already occurred in the vertical current sheet. At 
about t = 32t,, the velocity of the CMEs reaches the maximum 
speed, and the plasmas around the two current sheets are heated 
by reconnection (Figure 3(f)). Consequently, the prominences 
erupt as the main component of the three-part CMEs and are 
always accompanied by two-ribbon flares. 
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For case B, Figure 4(a) shows that the first catastrophe 
occurs at about t = 1574, and the length of the vertical current 
sheet below the flux rope is 0.02009 R,. Figure 4(b) shows 
that the second catastrophe occurs after t = 31r4, and then 
the prominences rapidly erupt into the interplanetary space. 
From 5t, to 20t,, the velocity of the rope has been within 
the range of 0.5kms^! to 1.0kms~! (Figure 4(c)). After 
a period of oscillation, the second catastrophe occurs, and 
then the prominences erupt and reach the maximum speed 
of 1090kms~! at 34t,4 (Figure 4(d)) Similarly, after the 
second catastrophe, magnetic energy of the system obviously 
decreases from 1.731 to 1.595 until the prominences obtained 
the maximum speed (Figure 4(e)), and the magnetic energy 
inside of the flux rope changes from 0.1338 to 0.05252 at 
the same time (Figure 4(f)). Furthermore, the corresponding 
magnetic field configurations are shown in Figure 5. At about 
t = 1414, the flux rope attaches to the surface of the Sun 
(Figure 5(a)), and then at about t = 15r4, the prominences 
form in the inner corona after the first catastrophe (Figure 5(b)). 
Before the second catastrophe, the prominences already levitate 
in the inner corona at about t = 31t,4 (Figure 5(c)). Due to 
the second catastrophe, the prominences break out; however, 
magnetic reconnection does not occur in the two current sheets 
at about t = 3274 (Figure 5(d)). Figure 5(e) shows that magnetic 
reconnection occurs in the vertical and transverse current sheets 
simultaneously. At about t = 3474, the velocity of the CMEs 
reaches the maximum speed, and the plasmas around the two 
current sheets are heated by magnetic reconnection (Figure 5(f)). 


4. CONCLUSION 


Previous work by Zhang et al. (2011) has concentrated on 
realizing that the physical features in the solar atmosphere, 
such as temperature and density have important influence on 
the evolution of the solar eruption event. In this paper, using a 
2.5-dimensional MHD simulation, we investigate the role 
played by the inner coronal null point in the formation and 
evolution of solar quiescent prominences. We find that mag- 
netic energy of the flux rope system is big when yy is small. 
The numerical results show that the bigger wy usually corre- 
sponds to the smaller catastrophic point, and the lower magnetic 
energy of the flux rope system or inside the flux rope. During 
the evolution of the flux rope system, while wy is greater, both 
the first and second catastrophe occur earlier, but the maximum 
velocity of CMEs with the erupted prominences is very closer. 
Thus, we confirm that the catastrophe occurs earlier and more 
easily when yy is bigger. However, yy has little effect on the 
maximum speed of eruption. 

Because there is no inner coronal null point in the dipolar 
field, vy is zero. We think that the prominences can be formed 
not in a dipolar magnetic field, but in a complex magnetic 
field with the inner coronal null point. Then, we argue that 
the topological structure with the inner coronal null point 
is a necessary condition for the formation and evolution of 
solar quiescent prominences. In a nutshell, it is easier for the 
prominences to form and to erupt as core parts of the CMEs in 
the structure with a greater Wy. 
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